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Overview 


The Projective Optimization Land use Information System (POLIS) is designed to 
provide land use, employment, housing and population forecasts for the subre- 
gional areas in the San Francisco Bay Area. Employment data is presented by 
five sectors --agriculture, manufacturing, transportation/FIRE, retail trade and 
services-- for each of the 107 zones. Additional outputs from the system include 
commuting flows by mode, auto and transit, shopping trips, while the model can 
be easily extended to estimate volume of retail sales. The model has been cali- 
brated with 1975 and 1980 data and has been used in the PROJECTIONS ’83 series 


of ABAG’s forecasts for the San Francisco Bay Area Economy (ABAG, 1983). 


POLIS allocates households and employment on the basis of several criteria. 
Residential choice is determined by the travel to work behavior and the availa- 
bility of housing, whereas the location of the retail sector is affected by the 
proximity to population centers and the inherent attractiveness of a zone for 
shopping. The locational patterns of the other sectors are influenced by the 
accessibility to labor, the existence and size of agglomeration economies and 
the structural characteristics of each sector. Land use controls are explicitly 
considered in the model and affect both the availability of housing and the loca- 


tion of economic activity. 


The proposed model is formulated as a nonlinear mathematical programming 
problem. Hence, the location of residences and housing and trip flows patterns 
are estimated in a single iteration and are consistent with the constraints of the 
model and each other. The objective function maximizes the locational surplus, 
utility, associated with the multimodal travel to work trips, the shopping trips 
and the agglomeration benefits accruing to employers. The linear constraints 


reflect the planning constraints, i.e. land use controls, exogenous location of 


part of the employment, and accounting identities --total employment to be 


allocated, origin/destination trip equations etc.--. 


POLIS completely replaces BEMOD and PLUM, the employment allocation and 
land use models for the Bay Area designed in the early ’70 ’s and the major 
source of all prior ABAG projections. POLIS differs from these models as well as 
most of the operational Lowry type models in several aspects. It considers in an 
integrated and internally consistent fashion the location of basic /nonbasic 
employment and housing. It reproduces the behavior at the micro level and 
analyzes the location of economic activity in terms of agglomeration economies 
and proximity to labour, rather than established patterns at the base year. 
Finally, it lends itself to the testing of issues such as land use controls, scarcily 
of housing and shifts in the sectoral composition of employment, issues which 
are expected to be the dominant policy concerns in the Bay Area during the 


Next Téen,years: 


As shown in Figure 1, POLIS is part of the the three tiered modeling system of 
ABAG. The Regional Information System (RIS), a dynamic input-output demo- 
graphic, economic and energy model estimates regionwide population and 
employment in 38 sectors (Brady and Yang, 1982). An econometric model, the 
County Employment Forecasting System (CEFS), allocates the jobs, aggreggated 
to 32 sectors, among the nine counties (Prastacos and Brady, 1983). Since the 
outputs of both the RIS and CEFS models act as control totals for POLIS, the 


latter provides a consistent set of subregional projections. 
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Introduction 


The evolution of urban areas is a complex phenomenon full of nonlinear and 
dynamic relationships. The latter arise from the specific nature of the process 
through which cities grow, mature and decline. In an urban environrnent there 
exists a multitude of actors--households, employers, developers, planners, etc.- 
-that respond to events from the outside world according to the objectives and 
constraints they face. The interaction of the different decision making units 


results into feedbacks that affect the way urban structure is formed. 


Urban modelling is the science that attempts to obtain a tractable mathernati- 
cal framework capable of reproducing these interactions. The field emerged in 
the early 1960’s out of the need to provide a tool which would assist the policy 
makers in their decision making and also to provide a greater insight into the 
theory of urban spatial structure formation. Over the past twenty years urban 
modelling has progressed remarkably and there are now several operational 


and/or theoretical urban models available. 


Lowry's urban activity model for Pittsburgh (Lowry, 1964) is probably the most 
well known of these. It has been the basis for most, if not all, of the operational 
models and has seen a considerable number of extensions (Garin, 1966; Goldner, 
1971; Wilson, 1974; Putman, 1977). It is a descriptive mode! attempting to cap- 
ture the interactions occurring between employment and population at a spa- 
tially disaggreggated level. Given the geographic distribution of basic employ- 
ment, population and service employment are allocated among zones using two 


spatial interaction functions. 


Lowry’s model and its derivatives have been criticized on several grounds. 


First, they lack a behavioral interpretation; the allocation process is not based 


on the concepts of utility or profit maximization but merely replicates aggre- 
gate trends. The functions used for allocating housing describe base year condi- 
tions and fail to address the decision making process of the individuals seeking 
residences. Second, they tend to disregard many of the economic activity/land 
use/transportation interactions. Each component of the system is considered 
in isolation of the rest and a rigid, often ad hoc, cause effect structure and top 
down information flow are assumed. Third, they are concerned only with the 
demand side of the land use system. Locators are permitted to locate wherever 
they desire since no constraints are imposed on the number or type of locators. 
Fourth, they do not provide a mechanism for handling planning constraints. 
The iterative solution of the model is not amenable to inequality constraints and 
despite the procedures suggested by Wilson (1971) the effect of land density con- 


straints can not be easily incorporated. 


The efforts of Mills (1967, 1972) and Ingram et. al. (1973) resulted in the creation 
of urban models which are more firmly rooted in the theory of urban economics. 
Mills follows Herbert and Stevens (1960) in proposing a linear programming for- 
mulation of the land use allocation problem. The difficulty with this approach is 
that it results in an efficient allocation of land uses which is not necessarily a 
realistic one. The allocation process is governed by a function representing 


aggreggate objectives (economic efficiency etc.) rather than individual behavior. 


The NBER model (Ingram et. al., 1972) is probably the most elaborate urban 
model designed to date. It simulates the housing market and consists of a set of 
submodels. The problems with this model are that its implementation requires 
data which are not usually available, that it does not attempt to simulate the 
location of the ernployment sectors and that only rudimentary attention is paid 


to the transportation sector. 


Ge 


More recently Coelho and Williams (1978) and Sharpe and Karlquist (1980) have 
presented a framework which overcomes some of the limitations associated with 
the Lowry and Mills models. Their models are formulated as mathematical pro- 
gramming problems and have been shown to generate Lowry-type results under 
certain assumptions. They are based on the concepts of random utility theory 
which provides a way of reproducing at the aggregate level individual behavior. 
Both of these models are theoretical and no attempt has been made to date to 


operationalize them. 


This study reports on the design and the empirical estimation of a planning 
model of this type for the San Francisco region. The proposed model, referred to 
as the Projective Optimization Land use Information System (POLIS), considers 
in a unified fashion the location of both basic/service employment and 
residences and can easily address planning constraints. Although the model is 
normative it is consistent with the behavior at the micro, individual level. Dis- 
tribution of trips and modal choice are integrated through logit functions, and 
the locational patterns of firms are influenced by the existence of agglomera- 


tion economies and proximity to labor force. 


The modelling effort was guided by three major objectives. First, to produce an 
operational model which avoids the pittfalls of PLUM -- the Bay Area implemen- 
tation of the Lowry model -- and encompasses some of the recent advances in 
the art of urban modelling, second to design a framework for which data are 
readily available or can be easily obtained and third to provide a tool which is 


amenable to the testing of planning policy issues. 


The set of output variables which can be obtained from POLIS include for each 
zone: 


1. Population 


102 


<. New housing units 

3. Employment, disaggreggated to 5 sectors 

4. Land use characteristics --vacant land, land consumed--. 

Oo. Number of work trips by mode and zone of residence and employment 

6. Number of shopping trips by zone of residence and location of retail 


establishments. 


Probabilistic Choice Framework for the Lowry Model 


The traditional Lowry model is characterized by a sequential decision making 
process; households first select employment location then look for a residence 
and finally select the location of other home based activities (shopping, recrea- 
tion, ete.). Ben Akiva et. al. (1978), Williams and Coelho(1978) and Wilson et. al. 
(1981) have proposed probabilistic models in which the residential choice is 
influenced not only by the location of employment but also from the shopping 


trips location. 


Within this framework the utility or surplus Uijmk of living at zone i, working 
at j, travelling to work by mode m and shopping at k is assumed to be 
represented by an additive separable function of the form: 

Ujjmx = Uij + Umi + Uxjijm (1) 
The first part Uj; defines the net utility of living at i and working at j and is 
affected by the characteristics of zone i as a residential choice, whereas Up); 
is the utility associated with mode m given that i and ] were chosen respectively 
as residential and employment locations. The last part Ux} ijm is the net surplus 
of shopping at location k given that zones i, ] and mode m were already selected. 
Modal choice on the shopping trips could be also introduced. Each of the net 
utilities can be defined as the difference between the benefits gained and the 


transportation or interaction costs incurred. All the net utilities are random 


= ies 


variables and therefore the theory of random utility (Domencich and McFadden, 


1975; Cochrane, 1975; Williams, 1977) can be readily applied. 


Assuming that Ux} ijm is independent of j and m, i.e. shopping choice is affected 


only by the residence location, we have: 


Une ap eter (2) 
yo ll Se W W 

Ce: Uhl ons e apmtrdes, (3) 

Ui; = bi SF wee = uve + ej (4) 


where uf is the trip end utility of shopping at k, ue, is the utility associated 
with mode m, see is the utility of living at 1, US is ameasure of the compara- 
tive advantage of shopping from i, us. is the comparative advantage (or disad- 
vantage) of commuting between i and j Cin is the perceived cost (or disutility) 
of travelling from i to j with mode m and cis the travel cost from i to shop- 


ping activities at k. Finally, es, enn and €j; are random variables represent- 


ing the variation in choice arising from unobserved factors. 


Because Uijmk is additive separable, the maximization of the utility Uijmx is 


consistent with the choice equation 


Pijmk = Pij * Pmjij ° Pklijm (5) 
where Pijmk is the probability of selecting zone i as a residential choice, j as an 
employment location, m as a mode of travel to work and k as a shopping loca- 
tion. Pjj is the probability of selecting i andj, Pmjij the conditional probability 
of choosing m given i and Jj, and PxJijm the conditional probability of shopping 
at k given that i, } and m have been already chosen. If the errors in (2), (3) and 
(4) are assumed to have a Weibull distribution and to be uncorrelated among the 
zones, then each of the components of (4) can be expressed as a nested multino- 


mial logit model. That is 


= De 


_exp [6° (ui = cf) ] , 
YexplA*(uf — ef.) (8) 


Pxlijm 72 
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exp[A* (ug — ch) ] 


 Dexp [Cus - ol] ‘iD 


where Be" ax erid ps are parameters related to the standard deviation of the 


error term in (2), (3) and (4). 


Omitting the Us term, i.e. disregarding any utility associated only with mode 
m and assuming that all the differences between modes are captured by the 


cost parameter cfm we have, 


W; WS exp (- BYus, a exp(—Ac#,) Wé exp (—B8c¥ 


Pajmk = EM, WH exp (Auf) “Yexp(—Ackn) | EWfexp (—#*oh) i. 
2 
where 
W; = exp (6 uy") (9) 
WS = exp(gs¥uS (10) 
WS = exp(f* ub (11) 


The functional form for the trip variables Tijm -number of work trips between i 
and j by mode m-, and 98;, -number of shopping trips between i and k-, is given 


in terms of total employment Ej and population H; as, 


Tym = Ey * Pi > Pmyjj = 


=8)¢ = oS OS 12 
J er Yue, +) Se oe o 
il m 


=fe3 5 


Wie exp (—BS ci) 


Sik = Hi > Pxyijm = Hi - Swe exp (—f8 cs, 
1 
k 


(13) 


The term us. is the mean surplus associated with commuting from i to j. 


Through (3) it can be defined to be 
1 
ue = - + InYexp(—Acijm) (14) 
m 


a result already obtained by Williams(1977) for the composite cost by consider- 
ing the generalized surplus integral of Hotelling (1938). Thus ule is the compo- 


site cost and can be replaced by the more commonly used symbol Cj; . 


Equations (12) and (13) define a model which has a similar structure to the 
results derived from the traditional Lowry model when a multimodal transport 
system is considered. The attractiveness weights of the Lowry model are 
shown to be nonlinear transformations of the trip end utilities. However, since 
it was assumed that the residential choice is optimized with respect both the 
work and shopping trips the residential attractiveness factor contains an addi- 
tional component ( Ws ) related to the comparative advantage of shopping from 


lL: 


The results obtained above, in addition to providing the Lowry model with a 
behavioral interpretation, permit the specification and estimation of a con- 
sistent measure for the benefits accruing from a specific policy. For example, 
the locational or consumer's surplus associated with a fixed distribution of { E; 


} can be defined as: 
Lov= >) Ej ujs (15) 
j 


where Uje the mean surplus associated with employment at j] is computed from 


(4) as 


mig Yee 


uje = ard exp[B¥(u* + us — or (16) 
1 
The location surplus can be then written as: 
LS = pi E. YS’ W: Ws (- Wo) 
aw yy pe iWsexp(—6" c;; fg) 
J 1 


and used for assessing the benefits arising from a change in the transport sys- 


tem (Neuberger, 1971). 


Mathematical Programming Framework 


As discussed earlier a major problem with the Lowry model and its derivatives is 
their inability to handle explicitly planning constraints. In this section we intro- 
duce a mathematical programming model which is shown to have a structure 
similar to that of the probabilistic choice model of the previous section and 


which can be easily extended to handle density constraints. 


Williams (1977) has shown that the total benefits associated with a spatial pat- 
tern of trips can be expressed in terms of the variables of a mathematical pro- 
gramming model. The objective function of this model is interpreted as the net 
utility of a group of individuals and its maximization reproduces the variability 
existing at the individual level. Coelho and Williams (1978) used this approach 
to analyze household location when it is influenced by work location and shop- 
ping trip behavior. The objective function in that case consists of the sum of two 
surplus functions, one related to the work trips { Tij } and one to the shopping 
trips { ij }. We will extend these results to accomodate a multimodal system 


(see also Boyce et. al., 1983). 


Consider a two sector, service and basic, representation of an urban economy. 
Assume that the residential choice is affected by the workplace location, the 


mode of travel to work and the shopping behavior. The mathematical 


at ae 


programming framework for a spatial interaction model is then given by: 


aga eA Go ; Sij ) = 


1 foe 1 
ak i BY arts ct eer emer at | mse } Tijml In Tim — 1) — ) Tijmei}m} + 
ijm ] 4m 1jm 
{ ‘i : [ in, ies s | (18) 
—_— —— om —_—_————— — es ‘Co 
B° | we up Say 
Ste 
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jm j 
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im 1 


where d, d’ and d”’ are scale parameters and EP is the basic employment in zone 
j. Constraint (19) ensures that the household population estimated from the 
work trip matrix is equal to that predicted from the shopping trips while (20) 
specifies that the total number of work trips to a zone is equal to the total 
employment (basic and service) in that zone appropriately scaled. The objective 
function consists of the sum of three surplus functions shown in (18) in brack- 
ets. The first two define the locational surplus attached to the residential 
choice when considering only work trips; the multiple dimensions of the deci- 
sion process associated with work trips, gives rise to two surplus functions, one 
for destination and one for modal choice. The last function depicts the loca- 
tional benefits arising when residential choice is also influenced by the shopping 
behavior. The parameters Be BE and A convert the utility associated with 
trip making into monetary units compatible with the transportation costs 


incurred. 


Puan 
To obtain the number of work trips Tijm consider the Langrangian: 
= Ae » 65 (4d Si; az » Tim) oF NII yy Sy ar d" EP rz bret, (22) 
1 J jm J i im 


where 6; and yj; are the dual variables of constraints (19) and (20). Since for 
finite a A and a the trip variables Tj; , Tijm and 5;j are always positive we have 


from the Kuhn-Tucker conditions; 


[ >) Tijm 
OT, neu | Be W, x nTijm) — Cijm — 55 — 77 (23) 
and 
i deta P| Sara gna (24) 
05ij pw we J 1 4) 


which can be rewritten in terms of the trip variables as: 


—r/ pW 
> Tim 
Tijm = ar? EXD =A Ci) EXP —A( 0; + 7) )) (25) 


Si) = WPexp[6*(d0; +d'7; — c#)] (26) 


Summation of (25) with respect m results to: 
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Replacing (28) into (25) we get: 


[ow 
Te i= |DexP(—Aciim) -r/ A+p" Sap Acie exp| egw (47) (29) 


Introducing a composite cost function cj; with form similar to that defined in 


«17s 


the probabilistic choice model we get: 
expC Aci) = ViexDCaAcinw (30) 
m 
and (29) can be rewritten after some manipulation as, 


exp (—Acijm) 
exp ( IN Shiaa) 
m 


Tijm = Wj’ exp 16° (: +75) exp (—B' cy) (31) 


eh 


where £' is defined to be equal to : 
p 4 BY +r 


Invoking the constraints for doubly constrained models, the balancing factors 


AM ; BY , AS and BP can be defined as: 


[ ni let 
ee | BYE; exp(-8 er (ae) 
J 
[ um™ al 
By = |» ANH; exp(—f ci) (33) 
Hf 
[ -1 
AP = » BeES exp 6 ae) (34) 
J 
[ $5) 
Be = y\ ASH; exp(—6 cj) (35) 
al 
and the trips variables expressed as: 
~ Exp =ACi) 
Tijm = Al” H; Bi’ Ej exp(—8 cy) een ade (36) 
ijm 
a J 
By > Ay Hi BP ES exp(—6"c3) (37) 


It can be easily demonstrated that the trip patterns (36) and (37) are of the 
same type with those obtained in the probabilistic model discussed earlier. 
Thus, the mathematical programming framework, although it has an aggregate 
nature, results in expressions at the micro level which are consistent with the 


theory of behavior at that level. Additionally, because of its optimization 
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nature, it can easily handle planning constraints. The objective function is con- 
cave separable subject to linear constraints and the model can be estimated by 


several available methods. 


The dual program associated with (18)-(21) is unconstrained and is given by 


min sg renee Pa = 


sei a =.) _exp(—Acim) 
= sr W, exp | fal Omit Cy, Sean 
B a iexp[— (6; + 7j + ci) ] TEEN 
m 
1 , ad 
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The objective function represents the total locational surplus accruing to loca- 
tors and suppliers of housing opportunities. The vector of shadow prices y is 
the consumer surplus of the individuals seeking a house, while 6 can be inter- 
preted as the location rent on housing. As Coelho and Williams{1978) point out 
estimation of the planning model must be carried out by solving the dual rather 
than the primal problem, since the number of variables in the dual is 


significantly smaller. 


Location of Employment 


In the model considered above only the retail sector’s location is determined 
endogenously. It was assumed that the location of the basic sectors is esimated 
outside of the model. With the exception to the work of Mackett (1980), this is a 
characteristic of all the land use models; basic employment location is 
predicted through a set of indexes and then used9ito “drive the 
residential/shopping choice model(s). Such an approach assumes that changes 


in housing location patterns will not influence the location of industries. 


aA19= 


Although this approach may be valid for the traditional heavy manufacturing 
sectors (oil,chemicals,transportation equipment) it may lead into erroneous 
results when considering the intraurban location of "footloose" industries which 
are highly dependent on access to labour markets. Electronics, Research and 
Development and several non retail services -- FIRE, part of Business and Per- 
sonal Services -- are typical examples of sectors where availability of qualified 
labor force and proximity to competitors are the major considerations for loca- 
tional decisions. It must be noted that employment in these sectors represents 
more than 50% of the non-retail employment af the Bay Area and is expected to 


grow at avery high rate. 


To integrate the location for all employment in the mathematical programming 
framework the factors influencing the locational decisions of the firms in the 
basic industries should be added to the criterion function. The most significant 
of these are: access to labor, access to markets, proximity to suppliers, 


agglomeration economies and land availability. 


The access to labor factor is already part of the model since the objective func- 
tion represents the employee’s consumer surplus with respect to their trip 
making and residential choices; therefore the accessibility of employment sup- 
ply to demand is maximized. Proximity to markets and supliers is of particular 
importance to retailers and industries with low value added per unit of weight 
ratio, but should not be a major consideration for firms producing for external 
markets or industries in the finance sector. For high tech type firms proximity 
to competitors/suppliers is desirable since it provides access to specialized 


labor. 


Agglomeration economies arise form the location of firms adjacent to each 


other. They include: access to sources of capital, labor market ecomonies, 


She 


communication and information economies, access to specialized business ser- 
vices, superior training facilities, attractive environment for innovation and 
other not easily quantifiable parameters. In the Bay Area they are highly 
significant and are the major cause for the rapid growth of the high technology 
industry in the Silicon Valley and the heavy concentration of financial and busi- 
ness services in the city of San Francisco. Unfortunately, they can not be 


estimated directly and must be defined in terms of surrogate variables. 


Agglomeration economies can be integrated into a land use model by separating 
them between those occurring at the macro level and those exhibited at the 
zonal level. The first ones consist of the causal relationships that link activities 
to each other. They represent comparative advantages and profitabilities aris- 
ing from the existing structure of production and consumption. The intersec- 
toral relationships shown by the Input-Output table is one manifestation of 
these. The propensity of certain industries to locate in specific areas is 


another. 


They can be incorporated in a mathematical programming model through the 
addition of equations in the constraint set linking employment in a particular 
sector to the size of the employment (lagged or current) in other related sec- 
tors and counties. These spatial sectoral equations can be defined similar to 


White and Hewings (1982) or Prastacos and Brady (1983) and take the form: 


Hie at (ame ny y, Ea s=ee) (39) 
ge 


where Q is the set of sectors with which strong economic relationships exist, t-1 
indicates a lagged value and the superscript * denotes regional variables. The 
choice of spatial level in (39) is governed by the availability of information. In 
order to capture long term trends the coefficients of (39) should be estimated 


from time series and not cross sectional data. SMSAs, counties or large cities 


Bes 


although they do not necessarily conform to economic boundaries, are 
appropriate levels for which long data series are usually available. In the pro- 
posed model for the San Francisco Area spatial sectoral equations similar to 


(39) are introduced for each county. 


The agglomeration economies at the zonal level arise from the inherent attrac- 
tiveness of certain zones. They are directly related to particular zonal charac- 
teristics and are sensitive to changes in these. Relative cost/availability of land 
and accessibility to other employment centers are probably the two major fac- 
tors that render some zones more attractive than others. Using microeconom- 
ics principles they can be interpreted as the marginal benefits accruing to loca- 


tors of employment opportunities. 


To incorporate them in the modelling framework assume that there exists a 
function iC) representing the agglomeration potential of every zone. Then, the 
component f(-)E could be introduced in (18) and the objective function be 
interpreted as the maximization of employer’s and employee’s surplus. The 
function f(-) can be an index or a linear/nonlinear transformation of zonal 


characteristics. 


The Land Use Information System for the Bay Area 


POLIS allocates employment and households on the basis of the criteria ela- 
borated earlier. Residential choice is related to the travel to work and shopping 
behavior, while retail activity is located with respect to residences and a factor 
representing the attractiveness of a zone for shopping. Manufacturing and 
other basic industries locational patterns are influenced by the accessibility to 
labor, the existence of agglomeration economies and the existing structure of 


production. 


OD 


The model is not dynamic, but it can be classified as quasi-dynamic since it 
simulates the change between two states. At each time period only the net 
increase in employment opportunities/households is allocated and relocation 
of base year jobs is handled by appropriately increasing the number of jobs to 
be distributed. Thus, the problem of creating an instant metropolis, for which 
the Lowry model has been criticized, is avoided. The total number of jobs and 
households to be allocated are given exogenously; hence, the results of the 


model are consistent with those derived from regional economic models. 


In describing the mathernatical framework of the system it is assumed that 
there are N zones and K economic sectors. The K sectors are partitioned into 
two sets, the basic sectors Kyas and the retail or local services sectors Kjg¢ . 
For variables with spatial and sectoral components subscripts identify zones 
and superscripts represent sectors. Base year values and control totals are 


denoted respectively with a superscript o and a bar. 


The complete model is written : 
max Z oon Shy Che AH; ) = 


23m 


1 
mas ijm i jm ijm 
1 Si n\a" apn 

pe a as ot WE oe Se ees Sib Gon ae Ve: (40) 

ke Koc k ij ii ne Koas 

subject to 

D Tim — ai (HP + AH;) = 0 (41) 
jm 
yin — aril Saiee a AE?) = =0 (42) 
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Eliade (He + AH;) = 0 (43) 
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Ss —h*(E#° + AEK) = (44) 
Sd "AE? < L, (45) 
n 

AH; < AH; < V; (46) 
), AE? SE a (47) 
J 

> AH; —-H.=0 (48) 
7 ae z pea Pe —y2=0 (49) 
AEP, < AEP < AEBy (50) 
Wee Saelad nO Ya (51) 


where AH; is the number of new households in i, AE;* is the number of addi- 
tional jobs\in Sectorm/)toybe Tocated at.j, L) is the area of land available for job 
location, si is the land availabie for housing at i and the rest of the variables 
take the definitions provided earlier. Py is the set of zones belonging in county 


Cc. 


The objective of the model maximizes the locational surplus associated with the 
multimodal travel to work trips, the shopping trips and the agglomeration 
benefits accruing to employers. Cost components are not included since the 
marginal cost of establishing a household or setting up a firm can not be 
estimated with any degree of accuracy. Equations (41) and (42) are the origin 
destination constraints for work trips and (43)-(44) fullfill the same function for 
the shopping trips. The coefficients es represent mean expenditure per house- 
hold of zone i for goods of sector k, hence, the flow variables si are interpreted 


as the volume of sales of goods of type k in zone j to residents of i. 


Constraints (45) and (46) are the land use constraints for employment and 


ay ee 


housing. Constraints (50) and the first part of (46) are introduced for some 
zones and sectors to account for exogenous location of a portion of the employ- 
ment (for example government) or other zoning considerations which are not 
captured in (45). Equations (47) and (48) define the total number of 
jobs/households to be allocated. Finally, linear constraints similar to (49) are 
introduced for each basic sector and county to reproduce established sectoral- 
spatial patterns and agglomeration economies at the macro level; they will be 


further discussed in the next section. 


The dual problem is written: 


EON? al Ai Opt. (eee Rmee. 5° Ten Rene 


1 pexp(= ACh) 
F a pL -8.(«; oy) exp ( =) Cijm) 
1 
Sas SN ema te ec.) | * RaiHy +L pp +t 
ij 


+ Sek ee ar early + y iVs Pye E Berle + bivoee 
ik 


such that 
bPo,; — dv, - f- gB+ Yes = (iP)e (53) 
q 
béo, + h*uk — dky; - & + DIkpd = 0 (54) 
q 
ayk; + ex ni ei oitey'0 (55) 


where Kj, Vj, nk, Lys, v}, pi, €, T, Po are the dual variables of constraints (41)- 


(49). 


There are two different ways of solving the proposed model. The first approach 
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is to solve the dual problem (52)-(55) using a nonlinear programming algorithm. 
The dual consists of a set of linear constraints and a concave function whose 
gradient can be readily estimated analytically at any feasible point, hence, any 
of the standard nonlinear programming algorithms (Conjugate Gradient, Convex 
Simplex, Quasi-Newton) can be employed. At optimality, the values of the dual 
variables uniquely define the primal ones and the optimal set of trip flows and 
population/employment allocation can be estimated as functions of the optimal 
dual variables (eq. (29), (26) and (42)). It must be stressed that, although the 
primal problem can be solved directly through a nonlinear programming algo- 
rithm --convex objective function subject to linear constraints--, the solution of 
the dual is preferable because the latter has a significantly smaller number of 


variables and constraints. 


The second approach exploits the structure of the primal. Since the nonlineari- 
ties arise from the trip variables and since the problem reduces to a trip distri- 
bution mode! when the housing and employment variables are fixed, a partition- 
ing or decomposition algorithm can be applied. Benders’ Partitioning Algorithm 
(BPA) (Benders, 1962; Geoffrion, 1972) is a natural candidate for such an applica- 
tion. It is an iterative algorithm which guarantees convergence to the optimal 
solution in a finite number of steps and which can provide a heuristic solution if 


computer time limitations force an early termination of the solution procedure. 


When applying BPA the model is partitioned into two problems, a master prob- 
lem which is a linear programming model estimating the location of housing 
(AH) and employment (AE) and a subproblem which is a combined 
distribution-modal split problem computing number of work and shopping trips. 
In the subproblem the trip origins and destinations take the values of the just 


computed AH and AE respectively. The solution algorithm consists of succes- 


2 - 


sively solving the master problem and the subproblem until convergence is 


reached. The dual of the distribution-modal split model gives rise to a con- 


straint 


related to trip costs and accessibilities which is then added to the mas- 


ter problem and the procedure is repeated until an acceptable level of conver- 


gence is achieved. 


The complete solution algorithm is outlined below and also shown in Figure 2. A 


more complete desription of the algorithm and its convergence properties can 


be found in Prastacos (1984). An index is added as a superscript on the left side 


of a variable to denote the iteration number. 


Step 0. 


Step 1. 


Select a convergence parameter €=O. Set the upper bound UB=+- and 
the lower bound LB=—-», Set the number of iterations N = 0. If a feasi- 


ble set of vectors ‘AH, !AE exists go to step 1; otherwise go to step 2. 


Increase N by one and compute the following combined distribution- 


modal split problem for fixed values of NAH and NAE. 


: d) lim : 
max W= — — eth Gn — 4) Bee tye me 
B™ ijm Wi ijm 
i i 


5 
— }} Tymeim = YF S $55 Cities 1) — Ske; (56) 
ijm k € Kice BK 5 W; ijk 


subject to the constraints (41) through (44). The two distribution prob- 
lems, one for Tijm and the other for Si, can be estimated indepen- 
dently of ech other. 

Denote the optimal value “W. Set LB=*r+%z,-—Nw ig 
LB= Nr +z, -—w. 

UB UB 


Test for convergence: lf UB = e€ the procedure has converged. 


Otherwise estimate the optimal dual variables k;, 9;, ni and us 


kOe 


Figure 2 
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associated with constraints (41) - (44) of the distribution problems and 


go to step 2. 


Step 2. Estimate the following linear programming problem. 


max le sAlZe tks ids (f7)? AED (57) 


in Ss Kas 


subject to: 


Lo = Y) (jai + np ey) AH; + 0) bP AEP + apn ARE” (58) 


i jn j,k € Kioc 
and constraints (45) through (51). 
C is a constant obtained from the dual of the distribution problems. Let 
N+1AH NHIAR and “48 denote respectively the optimal variables and 


value of this problem. Set UB = N+!F . Goto step 1. 


At the time this modelling effort was carried out we did not have access toa 
nonlinear programming computer package. Thus, the Bender’s Partitioning 
Algorithm was exclusively used for calibrating and estimating POLIS. A special 
computer prograrn was written to implement the algorithm and carry out the 
solution iterations. The linear programming master problem is estimated using 
XMP (Marsten, 1981), a hierarchically structured library of subroutines for 
linear programming, written and distributed by Professor Ray Marsten of the 
Department of Informations Systems at the University of Arizona. The com- 
bined distribution-modal split problem is computed using a standard entropy 


maximizing algorithm. 


BOGy 


Empirical Estimation of the Model 


The complete land use transportation model was applied to the San Francisco 
Bay Area (Figure 3) using aggregate data. The empirical estimation relied only 
on already available information, since lack of resources did not permit the 
undertaking of any major data collection effort, i.e. household surveys, special 
tabulations. Our experience with the implementation of the model indicates 
that models of this type can be succesfully estimated only when Census infor- 


mation are available. 


At the time this research took place the complete data sets for the 1980 Census 
were not yet released. As aresult, the model was calibrated using data from two 
different time periods 1975 and 1980. Information on the spatial distribution of 
jobs, housing and population were available for both 1975 - in house data base 
form a previous modelling effort - and 1980 - but trip table and trip cost data 
were available only for the year 1975. It is hoped that with the recent release of 
Urban Transportation Planning Package (UTPP) from the 1980 Census this 


discrepancy will be overcomed. 


The transportation data were obtained from the Metropolitan Transportation 
Commission (MTC), the transportation planning agency for the Bay Area. They 
included 1975 trip generation and attraction rates, detailed (440 zones) travel to 
work trip tables by mode, peak-hour level of service characteristics (travel 
time, distance) and values of time for the different modes. The trip table was 
‘synthetic’, in the sense that it was not produced from primary information but 
was obtained through MTC’s disaggregate travel demand modelling system (CSI, 
1980) and validated to replicate observed travel patterns at an aggregate level - 
30 zones- (Kollo et. al., 1981). Information on shopping trips was more scarce 


and consisted of an aggregate trip flows matrix and a detailed non-peak hour 
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level of service table. 


A special emphasis was placed in developing an uptodate accurate estimate of 
the amount of developed land and the supply of land. The Local Development 
Policy Survey Update, conducted in stages during 1981-82 by ABAG staff with the 
cooperation of local government agencies, revised the land use consumption 
data from the original 1976 ABAG Local Policy Survey. Development activity 
that occurred over the 1975-1980 period and post 1980 development potential, 
including both new development on vacant land and redevelopment, were 
assessed for each jurisdiction and merged with the earlier information. The 
resulting data base consisted of detailed land use profiles --developed land by 
type, vacant land and the potential for new jobs or residences -- for each census 


tract of the Bay Area. 


For the purpose of the model the nine counties of the Bay Area were divided into 
107 zones, each of the zones representing an aggregation of census tracts. Two 
modes, auto and transit, and four employment sectors were considered. The 
sectors and their SIC correspondence were: 

1. Manufacturing (SIC 20-39, 50-51) 

2. Transportation, FIRE (SIC 40-49, 60-67) 

3. Retail Trade (SIC 52-59) 

4. Services (SIC 70-89) 
This sectoral aggregation was necessitated in order to reduce the computational 
difficulties and most important to be able to produce forecasts at a level which 
is compatible with other studies of ABAG and MTC. lt is not the most appropri- 
ate one because it puts together sectors with different locational characteris- 
tics. For example, Transportation, a sector which is evenly distributed 


throughout the region, is lumped with FIRE, a sector which is highly concen- 


otp 


trated in the county of San Francisco. From the four sectors Manufacturing and 


Transportation/FIRE were assumed to be basic. 


In addition to these four sectors POLIS allocates employment in the Agriculture, 
Mining sector (SIC 01-14). However, since the employment in that sector 
represents a very small portion of the total regional employment and since its 
locational patterns are affected by natural resources, future employment levels 
are determined as a function of base year conditions and land availability rather 


than proximity to labor supply. 


The implementation of the model to the Bay Area consisted of three major 
tasks: estimation of the attractiveness weights and agglomeration economies 
function, specification of the spatial sectoral equations for the basic sectors -eq. 


(49)-, and finally calibration of the complete model. 


The attractiveness factors W; and W; represent the trip end utilities -eq.(9), 
(11)- and were defined as functions of different surrogate measures of attrac- 
tion. A major consideration on the definition of these factors was to estimate 
them in such a way so they can be easily reproduced in future time periods. In 
a forecasting environment these factors should not retain base year values but 
be updated every time period to reflect changing locational patterns. Addition- 
ally, if the model is to test alternative policies the attractiveness weights must 
be designed in such a way so that they are sensitive to the issues being tested. 
These are accomplished by defining them in terms of variables of the model or 


other variables which can be easily estimated outside the model. 


The attraction of a zone for residence was taken to be equal to: 


4 Hy 
See rum se (59) 


soo 


where H; and V; are respectively the already developed and vacant residential 
land and qd; is the ratio of median household income and median price of hous- 
ing at zone i. The q; ratios are normalized and can be interpreted as housing 
affordability indexes. They have the effect of rendering more attractive zones 
for which house prices are not high compared to incomes. Although q; can not 
be estimated for future time periods, we feel that being a ratio it is more time 


invariant than other housing characteristics like price or quality. 


For the retail and service sectors zonal attractiveness was defined as the pro- 
duct of the available land, the percentage of already developed land and relative 


accessibility to new households. That is, 


Y; EK 
WK =L, - Sf ee eee nee . ox 60 


where Yj is the total non residential developed land and Lj is the total avail- 
able land. The ratio id eZ Ej converts the percent of total developed land into 
land occupied by firms of sector k. The Hansen type accesibility index gi was 
introduced to indicate the propensity of local serving sectors to locate near new 


population centers. It is specified as: 
oe =u Alee wemne by cu) (61) 
i 


It can be shown that this definition is consistent with the rest of the model. The 
number of new households is lagged one time period to avoid estimation prob- 


lems. 


The absence of any information on the structure of the zonal agglomeration 
economies proved to be a major hurdle in the estimation process. There are 
infinite ways that they can be defined and no prior research work has be done in 
that field. Four different components were identified as the major components 


of an agglomeration function: size of employment at base year, rate of 


Sy \g 


employment growth ( sectoral and total) and proximity to other fast growing 
areas. Since the economies existing at the macro, county, level are captured by 
equation (49) the function f? must reflect relative intracounty differences 


appropriately normalized to account for the county and regional rate of growth. 


On the basis of these principles the function f? was specified as: 


ms Bret) SBP 7 Aegean Abi/ Boo sqy Gen 
*  —2> AES,/ AEB AR{,/ AE: AEP 


(62) 


where the subscript co denotes an employment variable for the county in which 
i belongs. All the AE variables are lagged one time period. Available land was 
not included in (62) since the agglomeration economies of a zone are more 


related to the existing rather than the future concentrations of employment. 


As discussed earlier the spatial sectoral equations are introduced to account for 
the agglomeration economies existing at the macro level and the structure of 
production and consumption. They describe the cumulative effect that these 
have on the location of industries rather than actually analyze the size or the 
composition of the economies per se. They were specified for each basic sector 


and county. For the Manufacturing sector they were defined as: 

Eeo =h ( Eéo,t-1 , AEs. , Heat) (63) 
and for Transportation, FIRE 

Beat = h ( ee , AEF .-1 , Bey , Hoe te eons (64) 
Subscripts co and * denote county and regional variables, t is the time period 


and the functions h and h’ are assumed to be linear. The employment E at any 


time period t is equal to the base year employment plus the change AE. 


The lagged variable is introduced to account for the historic concen‘ration of 


employment and the change in regional employment is a proxy for the rate of 


a FRESE 


growth occurring at the regional level. The other variables are included to iden- 
tify the cause and effect relationships between the different sectors. It was 
assumed that local services sectors are linked only with other activities in the 
same county while for basic sectors the linkages cross county boundaries. Thus, 
the employment variables of the former enter in the right hand side of (63) and 
(64) with their county values while those of the latter retain their regional 


values. 


Time series on county employment were obtained from the County Business 
Patterns reports for the year 1964-1980 and the equations (63) and (64) were 
estimated using ordinary least squares techniques. For alrnost all the counties 
and both sectors the fit was very good, Re greater than .90 and high significance 


values. 


The use of (63) and (64) for simulating the impact of the macro agglomeration 
economies may lead to erroneous results when producing long term forecasts. 
The coefficients of the regression equations can not be assumed to remain con- 
straint. This problem should be recognized when the model is used for forecast- 
ing activity levels for post 1985 years and the spatial sectoral equations should 
be partially relaxed. For the Manufacturing sectors a confidence interval of the 
form, 

(1-r)-h(-) s Ea < (1+r) ht) (85) 
could be defined. The variable r may take the values of .05 for year 1990 and .10 
for 1995 and 2000. An alternative procedure would be to specify confidence 


intervals for each coefficient in the regression; however, this would increase the 


number of constraints and further burden the computational requirements. 
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Calibration of the Model 


Calibration is the process of choosing values for the parameters of a model 
which produce the best correspondence between the model’s results and the 
observed situation at the base year. For the proposed model the variables to be 
calibrated are those which convert the utility associated with trip making into 
monetary units, namely chile BE, A and the exponents of the agglomeration 


function a”. 


The complex structure of POLIS --simultaneous residential/employment loca- 
tion and trip estimation --, and most important the lack of trip data for the year 
1980 did not permit the immediate application of any of the established calibra- 
tion techniques (Batty, 1976) for models of this type. Hence, an alternative cali- 
bration methodology was developed. The complete model was calibrated in 
three stages. The first two stages involved the calibration of the travel parame- 


ters, whereas at the last stage the a™ parameters were obtained. 


Initially, the model represented by equations (18)-(21) was calibrated on the 
1975 trip data. This is a model where only the location of the service sectors 
and residences are considered. It has been shown (Evans, 1971; Batty and 
Mackie, 1972) that for trip models with exponential cost functions the maximi- 
zation of the likelihood is equivalent to making the mean trip cost of Ae 
predicted trips equal to the actual mean trip cost. Thus, BY or actually 6 , 


BE and A can be estimated by solving: 
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There is a minor problem with this approach; the value of the left hand side of 


(66) is not known because the composite cost en is a function of the unknown 


robs 


parameter To facilitate the estimation of the three equations it can be 


assumed that cj = Coes 


An alternative approach avoiding this pitfall is to calibrate on the entropy of the 
transportation matrix. Erlander(1977, 1980), Chon{1982) and Boyce et. al.(1983) 
have shown that calibrating on the mean trip cost is equivalent to calibrating on 
the dispersion of the observed travel choices which is defined by the entropy of 


the system. In that case the three trip deterrence parameters can be estimated 


from: 
fa ie n im = fe In oi = oe (70) 


Rather than apply the commonly used Newton-Raphson method to solve (69), 
(70) and (71) simultaneously, a heuristic approach was followed. It was found 
that BE is not sensitive to changes in the work trip parameters, something 
probably attributed to the fact that the shopping trips matrix is more aggreg- 
gated than the one of the work trips. Thus, the shopping trips deterrence 
parameter is calibrated independently of the work trip parameters, while 
rc and A are estimated through a direct search procedure. The steps of the 


search procedure were the following: 


Step 1. Set number of iterations n=1 and assume an initial value of tt, 


ieee 


Step 2. Use "A to estimate Cij. Determine BM by calibrating on the distribution 
entropy Hp es d 
Step 3. Compute modal trip matrix Tijm and the modal entropy HY . Compare it 


with the observed entropy Hees 


. If they are equal the procedure ter- 
minates; otherwise go to step 4. 
Step 4. Estimate ae | through a three-point interpolation procedure similar to 


that first suggested by Hyman (1969) and later extended by I. Williams 


(1977). Increase the number of iterations n by one and go to step 2. 


As shown in Figure 4 this procedure was relatively efficient; it took twenty three 
iterations to reach an acceptable level of convergence. The calibration of ay 
(Step 2) was more rapid; by explicitly taking advantage of the fact that B<d 


(H. Williams, 1977) convergence was achieved in seven iterations on the average. 


Because the shopping trip data were not categorized by sector the same BE 
parameter was used for both Retail Trade and Services. An alternative pro- 
cedure would have been to calibrate on trip ends using the gradient of the likeli- 


hood function as the criterion of goodness-of-fit (Putman and Ducca, 1978). 


The second stage of the calibration procedure consisted of modifying the values 
for the trip parameters ol and BE so that the model reproduces the observed 
1980 distribution of housing and local services employment. Using 1975 as base 
year, the actual location of basic employment in 1980 and the calibrated modal 
split parameter A, the model was estimated for several sets of values of gw 
and BE. The set of trip distribution parameters calibrated on stage one were 
used to initialize the iteractive process. The goodness-of-fit criterion used was 
the R®* defined for housing as: 
SBUstgae tela a): 


SS (Hees ea H bs) (72) 
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and in a similar way for the service sectors. 


Finally in the last stage of the calibration the parameters on the agglomeration 
function were computed. Using the already obtained parameters rsh , A and 
gE and 1975 as base year, the model was estimated for several values of a”. 
Through a sequential search procedure the values of a” which maximized the 
correspondence, as measured by the R® criterion, between the predicted and 


actual locational patterns of the basic sectors were selected. 


Although the calibration methodology was heuristic, in the sense that it does 
not simultaneously estimate all the parameters, we believe that only marginally 
improved results would be obtained through a more vigorous approach. 
Through extensive testing it was determined that some of the parameters are 
mutually independent; for example both the a™ and the BE parameters were 
found to be relatively insensitive to small changes in the travel to work parame- 


ters. 


An interesting aspect of the calibration methodology is that it follows a bottom 
up approach. The transportation components of the model are calibrated before 
the locational part. This is unlike the approach used for most land use studies 
where a top down methodology results into calibrating first the parameters 


related to the locational decisions and then those describing travel behavior. 


The results of the calibrations are summarized in Table 1. The location of hous- 
ing and employment obtained by using the calibrated model to forecast 1980 
values are compared with the actual 1980 data. The goodness-of-fit of both the 
total values for housing/employment as well as the changes that occurred over 


the period 1975-1980 are measured. 


Overall the statistics depict a fairly good fit. The R® for total housing and for 


SAqie 


Table 1. Results of the Calibration 


Goodness of fit of predictions with actual 1980 data 


R 

Housing Units 

Total .90 

1975-1980 Change ve 
Total Employment 

All Sectors .89 

Manufacturing .90 

Transportation, FIRE 84 

Retail Trade 91 

services .86 
1975-1980 Empl. Change 

All Sectors si45 

Manufacturing 0 

Transportation, FIRE .64 

Retail Trade 82 

services .76 
Trips to Work (1975) 

Total .79 

Auto Bi 


Transit .69 
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each of the employment sectors are all between .85 and .91. The figures for 
total employment and housing .89 and .90 respectively indicate an almost per- 
fect fit. These high values also arise from the fact that the model allocates only 
the total regional change, which for the 1975-1980 period does not exceed the 


20% of the already in place housing/employment. 


The fit of the model when comparing the forecasts with the actual 1975-1980 
change was also acceptable. The R* for housing and employment drop to .74 
and .78. There is a wide variation in the fit of the different sectors. Retail Trade 
is the sector with best fit (.89) while Transportation, FIRE exhibits the worst 
(.64). The fit of the basic sectors is on the average less successfull then that of 
the retail/service sectors. This may be attributed to several factors: it is possi- 
ble that the zonal agglomeration functions were not defined or calibrated 
correctly, or that some of the factors influencing locational decisions of basic 
industries were ignored. Additionally, some errors are introduced from the way 
the basic sectors are defined; the sector of Transportation, FIRE includes 


employment groups which do not have the same locational characteristics. 


An interesting aspect revealed by the goodness of fit statistics is that the model 
is more successfull in predicting total employment rather than housing. Hous- 
ing which was not disaggregated either by ownership type or quality charac- 
teristics - price range or age and number of units in structure - has a fit which 
is not as good as the one for total employment. These results indicate that 
extensions in the model should be in the area of disaggregating housing by type 
and introducing supply equations linking supply and demand. The fit of total 
employment is also superior of that of most of the individual sectors, a sign that 


the model captures the aggregate locational patterns. 


When the results obtained at the 107 zone level are collapsed into the nine 


Aa 


counties and compared with the actual county changes the fit of the model 
improves noticeably. For all the sectors the R*® is higher than .85. The 
improvement is more dramatic for the basie sectors and is undoubtfully due to 
the existence of the spatial interaction constraints which relate at the county 


level the size of employment among to different sectors. 


The same improvement is also noticed when comparing the predicted and 
observed employment for twenty large industrial centers, like San Francisco, 
San Jose, Oakland-Berkeley, Walnut Creek-Concord, Santa Rosa etc. This 
reaffirms our belief that some of the problems are due to the way the zonal 


agglomeration function is specified. 


Only the fit of the travel to work trips is reported since there was not a detailed 
trip table for shopping trips. The fit of the work trips is not as good as expected. 
In aggregates transport studies R* in the range of .85 - .95 are not uncommon. 
However, it should be kept in mind that the model was calibrated to reproduce 
the 1980 distribution of housing/employment and not necessarily the trip table 
for 1975. The R®’s obtained after the first stage of the calibration process were 
close to .90. Finally, the generalized costs and values of time were taken from 


another study and might not be the appropriate ones for the proposed model. 


AAs, 


Use of the Model for Policy Simulations 


A major reason for undertaking the effort of building a new land use information 
system for the San Francisco Area was the need to provide a tool that would per- 
mit the testing of alternative public policies. Some of the policies which can be 
addressed by the model are: 


- What is the impact of changes in local policies regarding land 
development. 


- What is the impact of accelerated shifts in regional employment 
from Manufacturing to Research and Development, Finance and 
Services industries. 

- What is the impact of new investments in the transit system, 


like the proposed extension of BART, on the location of housing 
and employment. 


The first issue is related to the development policies of the different cities and 
counties of the Bay Area. A community’s development policies include general 
and specific plans and other programs to either encourage or discourage 
development activity in an area. Local zoning regulations for the type and den- 
sity of new developments, capital improvement schedules and building permit 


allocation are some of the methods used to manage the rate of growth. 


The different land use policies of local government can be used to define the 
supply of land available for accomodating future households and employment 
activities. Their impact can be simulated by the proposed model because land 
constraints are explicitly considered. For example, a policy imposing lower den- 
sities for new residential activities can be directly translated into number of 
acres available for development or potential housing units and its citywide and 
regiowide impacts on housing location be tested. More important, because the 
model addresses the issues of housing supply and employment location in a sys- 
tematic fashion, the impact of the policy restricting housing growth on the 


number of jobs attracted in the affected areas can be simulated. 
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The second issue, change in the sectoral composition of employment, can have 
profound effects in the character of the Bay Area. The new industries have loca- 
tional patterns and labor force requirements which are in several aspects 
different than those of traditional manufacturing. They are characterized by 
an increasing emphasis on decentralization and product specialization, reliance 
on the availability of a well educated labor force and diminishing requirements 
for large initial capital outlays. All of these might result into considerable shifts 
in population and employment. Employment in areas which traditionally have 
been considered to be ‘dormitory towns’ might suddenly swell because of the 
rapid increase in the number of self employed individuals and the formation of 


many small companies with specialized products, computer software for exam- 


ple. 


The emerging tendency in some industries of locating the office close to the 
homes of the employees might lead to a substantial divergence in the growth 
rates of the current CBDs and the peripheral areas; eventually, these will have 
important repercussions on the transport network utilization. Traffic volume on 
the CBD bound highways and transit systems might remain steady or even 
decline. In the peripheral areas --Central Contra Costa, Solano and Sonoma 
counties --, where auto is often the only available mode of transportation, the 


congestion in the highways might lead to chaotic situations. 


POLIS has the capability to simulate the impact that the changes in the compo- 
sition of regional employment will have at the local level and can assist in the 
evaluation of the various policy alternatives. Employment is disaggreggated in 
9 sectors and can be further disaggreggated should the need for a more detailed 
analysis arises. Transportation is integrated in the framework and trip flows 


matrixes for different growth scenarios can be computed and policy implica- 
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tions be derived. 


Finally the last issue, evaluation of investments in transportation infrastruc- 
ture, is the issue for which the majority of the land use models have tradition- 
ally been designed. Investments in the transportation infrastucture alter the 
accessibility of certain areas which in turn induce locational shifts. Addition- 
ally, as noted above, the change in the locational patterns might lead to an 
increased demand for transport in areas where the infrastructure is inadequate. 
It must be aknowledged that, because of the relatively high level of aggregation 
of the transport network, the model does not lend itself for the simulation of 
minor transportation investrnents, eg. highway interchanges, but only of major 
additions to the system. The proposed extensions of the BART systern and the 
provision of improved transit services in the high growth areas are examples of 


investments whose impacts could be tested by POLIS. 


Further Extensions 


There are clearly several areas where the model proposed in this study could be 
extended. The list of the potential extensions and further work can be aggre- 
gated into two major groups; the first group includes extensions to provide the 
model with an improved simulation of the housing markets and individual 
behavior, while the second group consists of procedures to calibrate the model 


using more efficient algorithms and more complete data sets. 


The goodness-of-fit statistics point out to the area where the major thrust of any 
extensions to POLIS must be focused. The simulation of the housing sector is 
not particularly successfull. Housing and household formation is assumed to be 
driven by the availability of land and the residential attractiveness factor. 


There is no disaggreggation by housing type or form of ownership and housing 
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prices are not considered. As a first step to relieve some of these problems, two 
types of housing, single family and multifamily, and two kinds of housing tenure, 
own and rent, must be introduced. The latter is particularly important because 
the characteristics of the households that own their residences is completely 


different than those that rent. 


As a second step, equations linking housing supply to housing demand must be 
designed. As formulated now, the model can not capture any relative changes in 
housing prices; it assumes that base year housing prices differentials will 
remain stable throughout the forecasting period. Housing supply equations 
similar to those used in the Chicago Area Transportation-Land Use Analysis Sys- 
tern (CATLAS) (Anas, 1983) or the NBER model can be specified. Land availability 
may continue to act as constraint to the provision of housing or be incorporated 
in the supply equations. In this case prices of new housing will be driven to 


infinity when the land supply has been exhausted. 


Household formation is another area where additional work is needed. Right 
now household sizes and number of workers per household are estimated exo- 
genously and are not tied to incomes. Introducing equations bringing together 
the size of the household, the number of workers, the type of tenure and the 
housing prices/rents will make the model sensitive to changes in the household 
formation rates and capable of answering questions with respect the demand for 


housing when the trend twords a smaller household size is reversed. 


The capability of the model to simulate the locational decisions of the different 
employrnent groups may be enhanced by redefining the zonal agglomeration 
economies function. Despite the introduction of the exponent in the f? func- 
tions, zonal agglomerations are almost linear and can not therefore capture 


some of the dynamic effects that take place through time. It may be worthwhile 
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to analyse the location of individual firms using random utility and replace the 
agglomeration function in the objective function by a component similar in 
structure to the one used for simulating the selection of residences by individu- 


als. 


There is also a need to disaggreggate economic activity into more uniform sec- 
tors. Although the four sectors categorization, used in the current implementa- 
tion of POLIS, is an improvement over the traditional basic/nonbasic separation 
of employment, important locational decisions are probably still hidden. Sec- 
tors with different growth rates and structural characteristics are lumped 
together. For example, the sector of Manufacturing represents Food Processing, 
a stagnant if not declining industry, and Electronics, a sector exhibiting a very 
rapid growth. Another factor to be considered when redefining the sectors is 
the impact of technology on the locational characteristics. The increasing 
importance of telecommunications and the change of the economy from one 
driven by the exchange of material goods to one where information is the main 


commodity will not affect all the industries in the same way. 


The adoption of the different extensions proposed above entails several techni- 
eal difficulties which must not be underemphasized. The introduction of non- 
linearities --housing and agglomeration economies -- must be limited to the 
objective function. Incorporating nonlinear constraints will render the model 
insolvable since there are no algorithms to handle large problems with non- 
linear objective function and constraints. Additionally, there is a tradeoff 
between the number of sectors and housing types considered and the difficulties 


encountered during calibration and solution time. 


The calibration procedure used in the current implementation of POLIS is a 


heuristic one and probably fails to capture some of the interdependencies 


= 4072 


existing among the various parameters. Theoretical integrity was sacrificed to 
the necessity of calibrating the model within a short time and a limited data set. 
A more complete calibration procedure based on the Newton-Raphson method 
must be designed. Also, additional work must be carried out on the tradeoffs 
between a simultaneous calibration of all the parameters and a calibration in 


stages. 


Finally the model must be calibrated again as soon as more complete data 
become available. The data used in the current calibration were from different 
time periods and different sources. The Bureau of the Census recently released 
the Urban Transportation Planning Package (UTPP) from the 1980 Census which 
contains important information on the travel patterns. It provides detailed 
commuting patterns by mode and census tract of residence and employment in 
addition to new socioeconomic characteristics. Use of this data set in leu of 
the 1975 trip patterns will result in a better simulation of the trip flows matrix 


and residence location. 
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Summary 


This study presented the structure and the empirical estimation of POLIS, the 
land use information and transportation model for the San Francisco region. 
The proposed model advances the state of art of land use modeling by overcom- 
ing some of the limitations associated with the Lowry framework. It is aggregate 
but it results in expressions which are consistent with individual behavior, it 
considers the location of basic/nonbasic sectors and residences in an 
integrated fashion, and can easily assess the impact of planning constraints and 


other development policies. 


In the last seven years the field of operational urban models has been in a stag- 
nant state. Disillusion with the early applications of the Lowry model has led 
most public agencies to abandon efforts to build large scale models. To our 
knowledge, POLIS is the only comprehensive land use - transportation model 
implemented for a metropolitan area of the United States in the last five years. 
The calibration of the model in the Bay Area points out that meaningful models 
can be calibrated even with limited data availability and limited resources. 
Finally, the successful use of POLIS for producing long range subregional fore- 
casts for the Bay Area is an indication that large models can be useful for plan- 


ning purposes if cast within the appropriate framework. 
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Guide to Notation 


Total eployment in zone j 

Total basic employment in zone i 

Regional employment in sector n to be allocated among the zones 
Employment in sector n in zone j 

Kmployment at the base year in sector n in zone j 

Housing units to be allocated among the zones 

Housing units (households) in zone i 

Total available land in zone } 

Locational surplus 

Probability of residing at i and working at j 


Probability of residing at i, working at j, commuting by mode m and 
shopping at location k 


Conditional probability of commuting by mode m, given that place of 
residence and work are respectively at zones 1 and ] 


Conditional probability of shopping at k, given that place of residence 
is at 1, place of work is at ] and mode m is used for commuting 


Number of shopping trips from zone i to service activities of sector k 
in zone j; or, volume of sales of goods of type k from establishments at 
zone j to residents of zone i 


Number of work trips from zone i to zone } 
Number of work trips from zone ito zone ] by mode m 
Net utility of living at i and working at j 


Net utility of living at zone i, working at zone j, travelling to work by 
moce m and shopping at k. 


Net utility of commuting by mode m given that place of residence and 
work are respectively at zones 1 and } 


Net utility of shopping at k given that place of residence is at i, place 
of work is at ] and mode m is used for the travel to work 


a 
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Vacant residential land in zone i 


Attractiveness factor for residential activity in zone i; represents trip 
end utility 


Attractiveness factor for service sector k in zone j; represents trip 
end utility 


Total non residential developed land in zone j 

Cost of traveling by mode m from i to j 

Composite travel cost between i and j 

Travel cost from i to shopping activities at location k 

Mean expenditure of households of zone i for goods of sector k 
Zonal agglomeration index for zone i and basic sector n 
Accessibility index for zone j and shopping activity k 


Housing affordability index for zone i; ratio of median household 
income and median housing price at zone i 


Mean surplus associated with employment in zone } 

Trip end utility of shopping at zone k 

Utility associated with commuting by mode m 

Utility of living at zone i 

Comparative advantage (disadvantage) of commuting from i to j 
Comparative advantage of shopping from zone i 

Number of new jobs of sector n in zone j 

Lower bound for the jobs of sector n to be allocated in zone j 
Upper bound for the jobs of sector n to be allocated in zone j 
New housing units in zone i 


Exponent for the agglomeration index of sector n 
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Deterrence parameter for the distribution of commuting trips; con- 
verts utility into monetary units compatible with the trip cost 


Deterrence parameter for the distribution of shopping trips; converts 
utility into monetary units compatible with the trip cost 


BY A 


Deterrence parameter equal to 


Dual variable of constraint (20) 
Dual variable of constraint (19) 
Dual variable of constraint (43) 
Dual variable of constraint (42) 


Dual variable of constraint (41) 


Deterrence parameter for the modal split of the journey to work trips; 
converts utility into monetary units compatible with the trip cost 


Dual variable of constraint (44) 
Dual variable of constraint (45) 
Dual variable of constraint (47) 
Dual variable of constraint (46) 
Dual variable of constraint (48) 


Dual variable of constraint (49) 
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APPENDIX A 


Input Data Needed for Using POLIS 


I. General 


. number of zones 

. humber of modes 

. number of employment sectors 

. housing units to be allocated 

. regional employment to be allocated 


II]. Calibration parameters 


. deterrence parameters for work trips (distribution and modal split) 
. deterrence parameter for shopping trips 
. agglomeration function exponent 


IJ. For every zone 


. base year population, housing units, employment (by sector) 
. base year land use conditions in acres 
- land use utilization by type (residential, basic + nonbasic) 
- vacant land by type (residential, basic + nonbasic) 
. densities 
- residential (housing units per acre) 
- employment by sector (employees per acre) 
. attractiveness factors 
- residential 
- nonbasic sectors 
. zonal agglomeration functions for basic sectors 
. workers per household, household size 


IV. Transportation 


. peak hour trip cost matrix for every mode (auto, transit) 
. non peak hour trip cost matrix for auto 
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Input Data for a Sample Zone 


Zone 20 Millbrae /Burlingame 


1985 Allocation 


I. Base year (1980) input data 


. Households: 24745 
. Total Employment: 44819 
Agriculture : 396 
Manufacturing : 5309 
Transp., FIRE : 23116 
Retail 5598 
Services : 10400 
. Developed residential land : 5068 acres 
. Developed basic + nonbasic land : 1439 acres 


II]. 1985 Densities and vacant land 


. Densities 
Residential : 4.92 units/acre 
Agriculture : 2.3 employees/acre 


Manufacturing : 30.7 
Transp., FIRE : 39.5 
Retail : 36.3 H 
Services : ge.d v 

. Vacant land (recycled land included) 
Residential : 27.6 acres 
Basic + nonbasic: 325.1 acres 


Il]. Attractiveness and agglomeration factors 


. Attractiveness factors (normalized) 
Residential : .037 
Retail Employment : .50¢ 
Service Employment : .810 
. Zonal Agglomeration factors 
Manufacturing : .053 
Transp., FIRE : .324 


IV. Other input data 
. Household size : 2.41 persons 


. Workers per household : 1.29 
. % Population in group quarters : 11.3% 
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